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Abstract: Non-local boundary conditions have been considered in theoretical high- 
energy physics with emphasis on one-loop quantum cosmology, one-loop conformal 
anomalies, Bose-Einstein condensation models and spectral branes. In the present pa¬ 
per, for the first time in the literature, the Wightman function, the vacuum expectation 
values of the held square and the energy-momentum tensor are investigated for a mas¬ 
sive scalar held satisfying non-local boundary conditions on a single and two parallel 
plates. The vacuum forces acting on the plates are evaluated. Interestingly, suitable 
choices of the kernel in the non-local boundary conditions lead to forces acting on the 
plates that can be repulsive for intermediate distances. It is then possible to obtain 
a locally stable equilibrium value of the interplate distance stabilized by the vacuum 
forces. 
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1. Introduction 


In recent years, non-local boundary conditions in quantum physics have been considered 
for at least three main purposes: 

(i) As part of the attempt of obtaining a consistent picture of one-loop quantum cosmol¬ 
ogy B.0 and in the course of investigating one-loop conformal anomalies for fermionic 
fields 0. 

(ii) Spectral boundary conditions for Laplace-type operators on a compact manifold 
with boundary are partly Dirichlet, partly (oblique) Neumann conditions, where the 
partitioning is provided by a pseudodifferential projection; they are of interest in string 
and brane theory ||, |5[ . 

(iii) As part of the investigation of bulk and surface states in Bose-Einstein condensa¬ 
tion models 0. 

In the latter case, following the work in Ref. [|J, it is useful to consider a simple 
example given by the Laplacian acting on scalar functions on the two-dimensional 
plane. More precisely, given the function q which is both Lebesgue summable and 
square-integrable on the real line, i.e. q G Li(R) fl ./^(R), one defines 0 

1 °° roc 

x) = E eilx/R J e~ ily/R q(y)dy. (i. 

l=— oo ’ x 



1 






















The function q r is, by construction, periodic with period 2itR , and tends to q as R 
tends to oo. On considering the region 

Br = {x,y : x 2 + y 2 < R 2 } , (1.2) 

one studies the Laplacian acting on square-integrable functions on Br, with non-local 
boundary conditions given by 


u -n 


+ f Qr( s ~ s')u(Rcos(s'/R), Rsm(s'/R))ds' = 0. 

dB R Jbb r 


(1.3) 


In polar coordinates, the resulting boundary-value problem reads as 

/ d 2 18 Id 2 


- 


y dr 2 r dr r 2 dip 2 


u = Eu, 


(1.4) 


du r 

R,<p) + Rj q R (R(<p-8))u(R,8)dO = 0. 


(1.5) 


For example, when the eigenvalue E is positive in Eq. 
function is 


(|1.4|) , the corresponding eigen- 


u hE (r,<p) = Ji(rVE)e lllp , (1.6) 

where Ji is the standard notation for the Bessel function of first kind of order I 6 Z. 
On denoting by q the Fourier transform of q, and inserting ( |1.6|) into the boundary con¬ 
dition (|E5D , one Ends an equation leading, implicitly, to the knowledge of the positive 
eigenvalues, i.e. 


VEJl(RVE) + Ji(Ry/E)q(l/R) 


0 . 


(1.7) 


The solutions of Eq. (|1.7| ) which decay rapidly away from the boundary are the surface 
states, whereas the solutions which remain non-vanishing are called bulk states [|6|. 

In the extension to gauge fields, non-local boundary conditions along the lines of 
Eq. (1.3) make it possible to improve the cllipticity properties of the boundary-value 
problem, by working with suitable symbols of the boundary operator, as shown in [:7|. 
The price to be paid, however, is that the gauge-field and ghost operators become 
pseudo-differential because the gauge-fixing functional is no longer a local functional of 
the gauge field [7j. 

On the other hand, in the field-theoretical analysis of macroscopic quantum effects 
such as the Casimir effect, an essential point is the relation between the mode-sum 
energy, evaluated as the sum of zero-point energies for each normal mode, and the vol¬ 
ume integral of the renormalized energy density. For scalar fields with general curvature 
coupling it has been shown that, in the discussion of this question for the Robin parallel 
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plates geometry it is necessary to include in the energy a surface term concentrated on 
the boundary ||. In subsequent work, by using variational methods, the first author 
of the present paper has derived an expression of the surface energy-momentum tensor 
for a scalar field with a general curvature coupling parameter in the general case of 
bulk and boundary geometries || . 

As a next step, we have been therefore led to consider the role of non-local bound¬ 
ary conditions in the course of studying the vacuum expectation value of the energy- 
momentum tensor as well as the Casimir effect itself. For this purpose, section [2] 
studies the Wightman function and Casimir densities for a single plate, while section 
[| is devoted to vacuum densities in the region between two parallel plates. Concluding 
remarks are made in section f|. 

2. Wightman function and Casimir densities for a single plate 

We consider a real scalar field <p(x) with general curvature coupling parameter £ satis¬ 
fying the field equation 


(V„V>‘ + m 2 + iR)<f = 0, 


( 2 . 1 ) 


where R is the scalar curvature for a (D + 1)-dimensional background spacetime, and 
V M is the covariant derivative operator. For special cases of minimally and conformally 
coupled scalars one has ( = 0 and ( = ( c = (D — l)/4 D, respectively. Our main 
interest in this paper will be the Wightman function, the vacuum expectation values 
(VEVs) of the field square and the energy-momentum tensor induced by a single and 
two parallel plates in Minkowski spacetime. For this problem the background spacetime 
is flat and in Eq. (|2.1|) we have R — 0. As a result the eigenmodes are independent of 
the curvature coupling parameter. However, the local properties of the vacuum such 
as energy density and vacuum stresses depend on this parameter. 

In this section we consider the properties of the vacuum for the geometry of a single 
plate. We will use rectangular coordinates x M = (t, x 1 = x, xy), where xy = ( x 2 ,..., x D ) 
denotes the coordinates parallel to the plate. We assume that the plate is located at 
x = 0 and the field obeys a non-local boundary condition similar to Eq. (|1.3| ), i.e. 



( 2 . 2 ) 


where n v is the inward-pointing normal to the boundary and the conditions on the 
function / will be specified below. For definiteness we consider the region x > 0 for 
which n v = 8”. It can be seen that, for this type of boundary condition, the scalar 
product (</?i, (f 2 )t for a given spatial hypersurface t = const, defined in the standard 
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way (see, for instance, [jpj) does not depend on the choice of hypersurface £. indeed, 
the corresponding difference for two hypersurfaces t = t\ and t = t 2 , from the field 
equation, by using the Stokes theorem, reads as 


(</W - (^1,^2)*! 


rt2 


= l 


't 1 


dt / dx||n y [<p*(f,0, X||)^(/7 2 (i,0, X||) - <p 2 (t, 0, xy)^^*^, 0, X||)] . (2.3) 


By virtue of the boundary condition (|2.2| ), the integral on the right-hand side vanishes. 
For the standard local Robin boundary conditions, the normal derivative of the held 
at a given point on the boundary is determined by the value of the held at the same 
point. The non-local boundary condition (|2.2|) states that the normal derivative at a 
given point depends on the values of the held at other points on the boundary. The 
properties of the boundary are coclihed by the function /. In some sense, the situation 
here is similar to that in electrodynamics for the spatial dispersion of the dielectric 
permittivity, when one considers the relation between the displacement and the electric 
held. In electrodynamics the spatial dispersion leads to the dependence of dielectric 
permittivity on the wave vector. Analogously, our non-local boundary conditions lead 
to the dependence of the coefficient F in the eigenfunctions on the wave vector ky (see 
below). 

As the hrst stage in the investigation of local quantum effects we consider the 
positive-frequency Wightman function. The VEVs of the held square and the energy- 
momentum tensor can be obtained from the Wightman function in the coincidence 
limit of the arguments with an additional renormalization procedure. Instead of the 
Wightman function we could take any other two-point function, but we choose the 
Wightman function because it also determines the response of particle detectors in a 
given state of motion. To evaluate the positive-frequency Wightman function we use 
the mode-sum formula 


(0 s |<p(^)<p(a; /M )|0s) = ( 2 - 4 ) 

k 

where |0g) is the vacuum state corresponding to the geometry of a single plate. For 
this geometry, the normalised eigenfunctions satisfying the boundary condition (|2.2| ) 
are given by 


gZk||X|| —iut 

<Ac( x M ) = , cos (kx + a), (2.5) 

V2 U 


where u 
[ 1:t- relation 


^ Jk 2 + k^ + m 2 , 0 < k < 00 , 


and the function a = a(k, k\\) is defined by 


2ia 

ik + F(k\\) ’ 


( 2 . 6 ) 
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with k\\ = I kill. In the last formula we have defined the Fourier transform 


F(h) = I dx,f (|x,|)e< k ii« 
(2tt) 


D — l 
2 


D — 3 

L 2 


d«« 2 f(u)Jp -3 


(2.7) 


where J w (x) is the Bessel function of first kind of order v. In the case F(k\\) > 0 there 
is also a purely imaginary eigenvalue k = iF(k\\) with the normalized eigenfunction 


vtV) = 


F(k ||) 


(27r) D - 1 a;( im ) 


D ik|| X|| —xF(fc||) 


( 2 . 8 ) 


arid = y/ k jj + m 2 — F 2 (A;|j), fc 2 > F 2 (fc||)—m 2 . These eigenfunctions correspond to 
the bound states of the quantum field. The occurrence of purely imaginary eigenvalues 
is proved by starting from the eigenfunctions, which depend on x according to A t e lkx + 
w pp constants Ai. From the boundary condition one finds eventually 

Ai(ik + F(k\\)) = A 2 (ik - F(k\\)). 

In addition to the solutions with both Ai ^ 0, this equation has indeed solutions 
Ai = 0 ,ik = F(k\\) and A 2 = 0 ,ik = —F(k\\). 

To escape the instability of the vacuum state, in the discuss ion belo w we will 
assume that the function ( |2.7| ) satisfies the condition F(k\\) < y^/c| + m 2 . For the 

convergence of the integral in (|2.7|) we need to have the behavior f(u ) = o(n 1_D//2 ) 
in the limit u —> oo and the behavior f(u) = o{u l ~ D ) in the limit u —> 0, by virtue 
of standard summability criteria at infinity and at the origin, respectively. For the 
discussion below of various asymptotic cases it is useful to have the behavior of the 
function F(k\\) for large and small values of the argument. For large values of k\\, by 
using the asymptotic expansion of the Bessel function for large values of the argument, 
from Q2.7|) we see that F(k\\) ~ o(k^ 


3-D 

2 


as fen 


oo. For small values of k\\ two subcases 


should be distinguished. For the case of functions f(u) ~ u —> oo, on using 

the expression for the Bessel function for small values of the argument, one finds to 
leading order 


Fn = 


2vrV 

h¥) 


duu u f(u), k\\ —>■ 0. 


(2.9) 


nh) 

For the functions f{u) ~ u — > oo, with D/2 — 1 < (3 < D — 1, by introducing 
in (|2 . 7| ) a new integration variable y = uk\\ and replacing the function f{y/k\\) by 
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Table 1: Examples of kernel function / in the boundary condition, with Fourier transform 
F\ _ 


f(x) 

F(v) 

.fo 

f (x 0 y k ( n 

(x 2 +X 2 )~ 2 v 

F r (Bjl-n) ^ y ) A n^y) 

/ 0 e -J?A A 2+x o 

(2 

f e~ v V* 2+ * 2 0 

op TT \?-l f K D,2-i(xoy/v 2 +y 2 ) 

U yj X 2 +X 2 

Z[Z7TXq) 2 fo ^2 +y 2)D/4-l/2 


its asymptotic expansion for large values of the argument, one finds F(k\\) ~ k^ +1 ~ D , 
—> 0. In the table below we present examples of the kernel function / in the boundary 
condition ( p.2|) with the corresponding Fourier transforms F. These functions depend 
on three parameters fo, xo, r/. 

The positive-frequency Wightman function is obtained by substituting the eigen¬ 
functions (|2.5|) and (|2.8|) into the mode-sum formula (|2.4| ). It can be presented in the 
form 


(0 5 |(p(x / ")(p(x /M )|0 5 ) = (0 M |<p(x^)<p(x' M )|0 M ) 

f rfkll iknfxii-x' ) f°° „ r , / /\ , 

+ / -— k|1 11 11 / dk -cos [k(x + x ‘) + 2a] 

J (2tt) J o u> 

r dk„ F(k l{ ) 

i k H(x||-x(|)-(*+ 3 : ')F(fc|,)+i(t'-t)« Clm) 

./ (2n) D ^ 


( 2 . 10 ) 


where |0 m) stands for the vacuum state in Minkowski spacetime without boundary, with 

the corresponding Wightman function, and 0 is the Heaviside step 
function. The last term on the right of this formula is the contribution of bound states. 
Note that, by using the definition for a, the cos function in the integrand can be also 
presented in the form 


n , o ^ ^-F\k l{ ) 1 2kF(k ll ) . , 

cos (ky + 2a) = — T _ i9n \ cos ky — yy —— ^ on N sin ky. 


( 2 . 11 ) 


k 2 + F 2 (k,) “ k 2 + F 2 (k t ) 

The integrals in (|2.10 ) over the angular part of the vector ky can be explicitly evaluated 
with the help of the formula already used in (|2.7|) . In the coincidence limit from formula 
(|2.10| ) one hnds the VEV of the held square 


( 0 s|<p 2 | 0 s) = + 


Sd -i 

( 2 ^r 


D 


dk\\ k\ 


D—2 


dk 


cos (2 kx + 2a) 


UJ 


■ . S °~ 1 1 [ dh\k^- 2 ^^d(F(k\\))e~ 2xF{ - k \\\ 

(2tt) d - 1 J 11 II u/™) V v 11 


( 2 . 12 ) 
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where Sd-i = 2 / rr l ' D ~ 1 ^ 2 /T((D — l)/2) is the surface area of the unit sphere in D- 
dimensional space. The last two terms on the right of this formula are induced by 
the plate. For points away from the boundary {x 7 ^ 0) these terms are finite and the 
divergences in the VEV of the held square are contained in the first term only. Hence, 
here the renormalization procedure is the same as that for Minkowski spacetime without 
boundary. 

To obtain an alternative form for the VEV of the held square, we write the cos 
function in the second term on the right of (| 2 . 12 |) in terms of the exponential functions 
and rotate the integration contour in the complex fc-plane by the angle 7 t/ 2 for the 
term with e 2lkx and by the angle —7 t/2 for the term with e~ 2lkx . We assume that the 
points ±i^Jk 2 + m 2 and the poles ±iF(k\\) in the case F(k\\) > 0 are bypassed on the 
right by semicircles with small radii. In such a way the following relation is obtained: 



dk 


cos ( 2 kx + 2a) 
a; 


-.—2 ux 


du 


u + F(k\\) 


k„+m 2 


V 


-2nF(k\\) 


u 2 — k 2 — m 2 u F(k\\) 

- 2xF (h)e{F{k ll )) 


^Jk 2 + m 2 — F 2 {k \\)' 


(2.13) 


Now we see that the second term on the right of this formula cancels out the third term 
on the right-hand side of formula (| 2 . 12 |) and we obtain 

Sd-i 


(p )sub = (0s|<£ 10g) — (Om\p |0m) = 


D 


n —2ux 


X 


du 


( 2*0 

u + F(k\\) 


dk\\ k. 


D—2 


v k l +m2 


k 2 


m 


2 u — F{k\\) 


(2.14) 


The VEVs of the held square in the cases of Dirichlet and Neumann boundary con¬ 
ditions are obtained from the general formula ([2.14|) in the limits F{k\\) —> 00 and 
F{k\\) —> 0, respectively. In these cases the integrals are evaluated by introducing a 
new integration variable v = ^ ju 2 — k 2 — m 2 and passing to polar coordinates in the 
(fc||, u)-plane. This simple calculation leads to the result 

(. m/x ) £ i li 

{{p )sub) Dirichlet = - ({p )sub) Neumann = d +1 K( D -i)/2{2mx ), (2.15) 

2 U 7T 2 


where K v {x) is the modihed Bessel function of second kind. The VEV of the held square 
induced by a single plate on the background spacetime x E, with an internal 

space E and local Robin boundary condition, is investigated in [lTJ as a limiting case 
of the braneworld geometry. 
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Having the Wightman function (|2.10|) and the VEV of the field square we can 
evaluate the VEV of the energy-momentum tensor by making use of the formula 

(OsliyOs) = Uni(0 s |a^(i“)9^(i'")|0 s ) + (c - Ij g„Ad“ - 0A (0sl»j' 2 |0s). 

(2.16) 

With the help of formulae ( |2.10| ) and ( |2.14| ) for the boundary induced part the following 
result is obtained (no summation over u): 


ecu = <0sK|0 s ) - (o M |rr|o M > = 


' M 
r»oo 


Sd-i% 

(2tt) d 


dk\\ k 


D—2 


X 


du 


A v (u,k\\)e 2ux u + F(k\\) 

l^/kj+m 2 \j u2 ~ kf\ ~ 777,2 U ~ F{k\\) 

where we have introduced the notations 


(2.17) 


Aq(u, k\\) = k 2 + m 2 — 4 (u 2 , Ai(u, fc||) = 0, (2.18) 

k 2 

A„(u, k l{ ) = (1 - 4C)u 2 - u — 2,..., D. (2.19) 

For the evaluation of the T" components, with v = 2,3,..., we note that, from the 
problem symmetry, one has T 2 = = ... and hence 




1 V 


P-1)' 


For the sum in the last formula, the integrand contains the factor 


D 

J]r;^-q + « 2 (i-4C)(D-i). 

12=2 


As we see, the vacuum stress in the direction perpendicular to the plate, (T)) sub , 
vanishes. This result could be also directly obtained from the continuity equation 
dv(T £)sub = 0 for the boundary induced VEV of the energy-momentum tensor. It can 
be checked that the tensor (|2.17|) is traceless for a conformally coupled massless scalar 
field. In the case of the local boundary condition of Robin type the function F(k\\) is a 
constant and the expression for the corresponding energy-m omentum tens or is further 
simplified by introducing a new integration variable v = yJu 2 — k 2 — m 2 and passing 
to polar coordinates in the (fey, u)-plane as before (2.15). In particular, for the massless 
case we obtain the result given in [§]. Note that in this case the VEV of the energy- 
momentum tensor vanishes for a conformally coupled scalar field. For the non-local 


















boundary condition, in general, this is not the case. Another important difference is 
that, unlike the local case, for non-local boundary conditions, in general, (Tg) sub 7 ^ 
(T>) su b (no summation over // = 2 For Dirichlet and Neumann boundary 
conditions the vacuum energy-momentum tensor is further simplified by a method 
similar to that used for the VEV of the field square, and one finds (no summation over 
V) 


«iy>. *) 

(m/x) 
2 D n 1 ^ 1 


Dirichlet 
D+l 


' Neumann 


= - «r;),ub) 

2 — Cc)Kd± 1 (2mx) + (4£ — l)mxKp~i (2mx) 


( 2 . 20 ) 


H = 0, 2,..., D. I 11 the Dirichlet case the corresponding energy density is negative 
everywhere for both minimally and conformally coupled scalars. The vacuum energy- 
momentum tensor for a plate on the background spacetime R ( ' D,1 ' > x E, with an internal 
space £ and local Robin boundary condition is investigated in Ql6 |. 

Now let us consider the limiting cases of the VEVs (|2.14|) , (|2 .1 71 ) . These VEVs 
diverge on the boundary. Surface divergences are well-known in quantum field theory 
with boundaries and are investigated for various types of boundary geometries and 
boundary conditions (see, for instance, IBB). They result from the idealization of the 
boundaries as perfectly smooth surfaces which are perfect reflectors at all frequencies. 
It seems plausible that such effects as surface roughness, or the microstructure of the 
boundary on small scales can introduce a physical cutoff needed to produce finite values 
of surface quantities. An alternative mechanism for introducing a cutoff which removes 
singular behaviour on boundaries is to allow the position of the boundary to undergo 
quantum fluctuations Jl3|. Such fluctuations smear out the contribution of the high- 


frequency modes without the need to introduce an explicit high-frequency cutoff. Note 
that in this paper we consider boundary induced vacuum densities which are finite 
away from the boundary. We expect that similar results would be obtained in the 
model where instead of externally imposed boundary condition the fluctuating field is 
coupled to a smooth background potential that implements the boundary condition in 
a certain limit [[14 1. I 11 the problem under consideration, for the points near the plate, 
the main contribution results from large values of u. To leading order we can omit 
F(k\\) in the integrands and the VEVs behave as those for the Neumann boundary 
conditions, i.e. 


r i -d 

<^ 2 >sub « ——D+yr 

(47t) 2 


D — 1 



2C(Cc-C) r 

( 47 r) 2 x D+1 


D+l 


( 2 . 21 ) 


For large distances from the plate the main contribution to the VEVs comes from 
small values of u. It can be seen that, for a massless field or for a massive field with 
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/(«) ~ u /3 , D /2 — 1 < /3 < D — 1, at large distances the VEVs behave as those for 
Dirichlet boundary condition. For the case of massive held with f{u) ~ o(u 1 D ), to 
leading order one has 

{y )sub ~ ^ _ m ((x )sub) Dirichlet , X —> oo, (2.22) 

and a similar relation for the VEV of the energy-momentum tensor. Here F 0 is dehned 
by relation (pT9|) . 

We have done numerical calculations for the components of the vacuum energy- 
momentum tensor by taking, for example, the kernel function 

f(x) = foe~ vx . (2.23) 

The corresponding Fourier transform F(k\\) is obtained from the second line of the table 
given earlier in the limit Xo —> 0, i.e. 

F (k,) = 7 -< 2 ’ 24 > 

f 1 + 

with the notation F\ = 2 D ~ 1 7T^~ 1 T(D/2)f 0 /r] D . In hgure 1 we have presented the 
vacuum energy density (full curves) and Tf-stress (dashed curves) as functions on 
rjx for various values of the parameter F\ (numbers near the curves) for a minimally 
coupled D = 3 scalar held. The energy density is positive for the region near the plate 
and is negative at large distances from the plate, having the negative minimum for 
some intermediate distance. The corresponding curves for the energy density of the 
conformally coupled scalar held are presented in hgure 2. In this case the Tf-stress is 
related to the energy density by the traceless condition: (T|) sub = —(TQ) sub /2. Note 
that, for the case of a conformally coupled massless scalar held with local boundary 
condition, the corresponding vacuum energy-momentum tensor vanishes. 
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1 1.5 2 2.5 3 3.5 4 


Figure 1: Single-plate induced vacuum densities (TQ) su b/?? D+1 (full curves) and 
(Ti) sub/v D+1 (dashed curves) as functions of ijx for a minimally coupled massless scalar 
in D = 3 in the case of kernel function ( [2.23 ). The numbers near the curves are the values of 
the parameter F\ defined in the paragraph after formula ( 2.24 ). 
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Figure 2: Single-plate induced vacuum energy density (Tq) su b/ r] D+1 as a function of rjx for 
a conformally coupled massless scalar in D = 3 in the case of kernel function ( |2.23 ). The 
numbers near the curves are the values of the parameter F\ defined in the paragraph after 
formula ( 2.24 ). The vacuum stress (T|) su b is expressed through the energy density by the 
traceless condition. 
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3. Vacuum densities in the region between two parallel plates 

3.1 Eigenfunctions 

In this section we investigate the VEVs for the geometry of two parallel plates. We 
assume that the plates are located at x — a,j, j — 1, 2 and the field obeys non-local 
boundary conditions 

n v {j) d u cp(x^) + J (fxj| fj (|x|| - xj||)^(a; /M ) = 0, x = a j: (3.1) 

where njV is the inward-pointing unit normal to the boundary at x = a r Here we 
consider the general case when the kernel functions fj determining the properties of 
the boundaries are different for the plates. The VEVs in the regions x < cq and 
x > <22 coincide with the corresponding quantities for a single plate geometry and are 
investigated in the previous section. Here we consider the region between the plates 
with njV = (—1 In this region the corresponding eigenfunctions are presented 

in two equivalent forms corresponding to j = 1, 2, i.e. 

= /3(k, /c||)e* k|iX||-4a ^ cos [k\x — a,j\ + ctj], 

where oq = aij(k, fey) is defined by the relation 

2iaj(k,k\\) _ ~ ( — J|) 

~ ik + i-iy-^Fj^y 

with the definition of Fourier transform Fj(k\\) similar to Eq. ( |2.7p , i.e. 

Fj(k ,) = 

duu~^~ fj(u)Jp- 3 (uk\\). (3.4) 

The corresponding eigenvalues z = ka , a = a 2 — oq, are solutions of the following 
transcendental equation: 

( z 2 — cic 2 ) sin z + (ci + c 2 ) z cos z — 0, (3.5) 

where the coefficients Cj are dehned by the relations 

Cj = l.-iy-'aFjk,). (3.6) 


( 2 < 


D — l 
2 


D — 3 

k, 


(3.2) 


(3.3) 
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Equation ( |3.5|) is obtained by bearing in mind that the ^-dependence in the eigenfunc¬ 
tions has the form Aie lkx + A 2 e~ lkx . The boundary conditions ( [3.1|) lead therefore to 
the linear homogeneous system 


(Fi + ik)e ikai A 1 + (F 1 - ik)e~ ikai A 2 = 0, (3.7) 

(F 2 + ik)e ika2 A l + (F 2 - ik)e~ ika2 A 2 = 0. (3.8) 


Non-trivial solutions for A 1 and A 2 exist if and only if Eq. 


holds. The expression 


for the coefficient (3(k, /q) in (|3.2| ) is obtained from the normalization condition: 


/3 2 (ka,k\\) = (2nf 1 au 


sin (ka) 

1 H--—- cos (ka + 2aj) 

rvCL 


(3.9) 


The eigenvalue condition (|3.5|) has an inhnite set of real zeros which we will denote by 
z = X n , n — 1,2,.... In addition, depending on the values of the coefficients Cj, this 
equation has two or four complex conjugate purely imaginary zeros (see, for instance, 
0) ±iyi, Vi > o. 


3.2 Wightman function 

Substituting the eigenfunctions ( |3.2| ) into the corresponding mode-sum formula, for the 
positive-frequency Wightman function in the region between two plates one finds 

(0|<p(rr /i )<p(;r /A1 )|0) = / dk||e lk||( ' X|l_x ii^ —— cos(zXj + aq) cos(zx'j + aj), 


Z=\n,iyi 


P 2 (z,h) 


(3.10) 

where |0) is the vacuum state in the region between the plates and we use the notations 


\x — a,- 


\x — a, 


Xj = 


-, Xj = 


a J a 

The summation over the eigenvalues A n , iyi can be done by using the formula 

h(z) 1 h( 0) 1 f0 ° 


(3.11) 


E 


' 1 + cos(x: + 2ai) sin z/z 

Z=\n,iyi 


2 1 - Cl 1 - C 1 


7T 


dzh(z) 


dt 


n Jo 

2 c,- 


h(te m / 2 ) — h(te m ^ 2 ) 

(t~ci)(t-c 2 ) 2 t _ i 

(t+Cl)(t+C2) 


[^( Ci e-/ 2 ) + ^( c .e-W 2 )], (3.i2) 


where 


9j(z) = (- 2 + c i) h ^)‘ 


(3.13) 
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Formula (|3.12|) is derived in [§ as a special case of the generalized Abel-Plana summa¬ 


tion formula |T5|. As a function h (with first-order poles at z = ±icj) we take 

e iw(t'-t) . - 

h(z) =- cos(zXj + a.j) cos [zx) + aq), u ; = J z 2 /cl 2 + k) + m 2 , (3-14) 

CLUJ V 


where the function a 3 = aj(z/a , k\\) is defined by formula (|3.3| ). By using the relation 
(from Eq. (|2 ■ 1 1|) ) 


~ 2 — c 2 2 zc ■ 

cos (y + 2aj) = J cos y -sin y, 

Z 2 + C 2 z 2 + 


(3.15) 


it can be seen that 


9 j{cje m/2 ) + gj(c je m/2 ) = 2c) 




aui 


(im) 


g-c 3 -(xj+x'-)_ 


(3.16) 


Moreover, by making use of the definition for aq we see that e 2ia A 0,fe n) = — 1, and hence 
cos(2«j(0, fc||)) = —1. This implies in turn that h(0) = 0. The resulting Wightman 


function from (|3.10|) is found to be 

<0|V2(£>(^)|0) = (0 S |^)^(^)|0 5 ),- + 

/»oo 

dt 


dkne 


ik n( x ir x Ii) 


x 


3 ' ( 2 tt) d j 
cosh (far,- + dij) cosh(hx) + ctj) 


(t-g)(t-c2) 2t _ i 
P+ei)(t+C2) C 


X 


J aJk 2 +m 2 
cosh (t — t')yjt 2 /a 2 — k) — m 


t 2 — kla 2 — m 2 a 2 


(3.17) 


where the function dtj = &j(t, /cm) is defined by the relation 


2 a ■ ^ O 

e 3 = J 


t + Co 


(3.18) 


and (0s|</j(x At )</?(x ,/1 ) |0s)j is the Wightman function for a single plate located at x — aj. 
On taking the coincidence limit, for the VEV of the field square we obtain the formula 


(0|<p 2 |0> = (0slv? 2 |0s>,+ 


4Sp-i 
(2tt 


i D 


dknk. 


D—2 


dt 


COsh 2 {tXj + OLj) 


a \J k l+ m ‘ 2 ^Jt 2 — k 2 a 2 — m 2 a 2 


x 


(t-Ci)(t-C2) c2t _ 1 
_ (t + Ci)(t + C 2 ) 


n -1 


(3.19) 
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where (Os|</? 2 |Og)j is the corresponding VEV for the geometry of a single plate at x — a,j 
and is investigated in the previous section. The surface divergences on the plate at 
x = cij are contained in this term. The second term on the right of formula (|3.19|) is 
finite at x = dj and is induced by the second plate located at x = a n , jT = 1, 2, j i ^ j. 
This term diverges at x = a n . The corresponding divergence is the same as that for 
the geometry of a single plate located at x = a n . 


3.3 VEV of the energy-momentum tensor and vacuum forces 

The vacuum expectation value of the energy-momentum tensor is evaluated by formula 
(|2.16| ) with the vacuum state |0). By taking into account formulae ( |3 • 1 7| ) , (|3.19|) , for 
the region between the plates one finds (no summation over u) 


2 S' 

(o|r;|o) = (o 5 |r;|Og) J + ^^i 


dh\k i 


D—2 


k\\ 1 x)dt 


'a^Jk^+m 2 ^ jt 2 — k 2 a 2 — m 2 a 2 


x 


(t - Ci) (t - c 2 ) ' 

(t + Ci)(t + C2) 


,2 1 


(3.20) 


where (OslT^IOg)^ is the vacuum energy-momentum tensor for the geometry of a single 
plate located at x — aj, and the second term on the right is the part of the energy- 
momentum tensor induced by the presence of the second plate. In formula ( |3.20|) we 
have defined 


k\\,x) = fcjj + m 2 — t 2 /a 2 + ( k jj + m 2 — 4 (t 2 /a 2 ) cosh (2 txj + 2 aj ), 
fji(t,k\\,x) = t 2 / a 2 , 


fjv{t,k\\,x) 


k 2 

W^T) 


k 2 

(D^T) 


+ (4C - 1) 



cosh (2 txj + 2 aj), 


(3.21) 

(3.22) 

(3.23) 


with v — 2, 3,..., D. It can be easily checked that the vacuum energy-momentum 
tensor is traceless for a conformally coupled massless scalar field. As we could expect 
from the problem symmetry, the vacuum stresses in the directions parallel to the plates 
are isotropic. Note that the T^-component of the energy-momentum tensor is uniform. 
This also follows from the continuity equation for the VEV of the energy-momentum 
tensor. In particular, the T^-stress is finite on the plates and determines the vacuum 
force acting per unit surface of the plate: 


o q roc 

P = -< 0|I?1 0) = wi 

'{t-F 1 (u)){t + F 2 (u)) 
X .(t + («))(;£-F 2 (u)) 


ro o 

duu D ~ 2 j 

J Vu‘ 2 -\-m ' 2 



t 2 dt 

yjt 2 — u 2 — m 2 


(3.24) 
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When the functions Fj(u) = const, this formula can be simplified by introducing a new 
integration variable v = \/t 2 — u 2 — m 2 and passing to polar coordinates in the (u, v ) 
plane. After integrating over the angular part one hnds the formula 


P = 


2 1 -d 7T -d/2 

WM 


dt 


t 2 ( t 2 - m 2 ) 2 

(t-F 1 ){t+F 2 ) 2at 

o t+F 1 )(t~F 2 r 


-1 


(3.25) 


For a massless scalar field this result coincides with the formula obtained in ||. The 
VEVs of the field square and the energy-momentum tensor for the geometry of two 
parallel plates on the background spacetime x E with an internal space E and 

local Robin boundary conditions are investigated in [11], |1|. The vacuum forces for 
Dirichlet and Neumann boundary conditions are the same and are obtained from (|3.25|) 
as special cases with Fj = oo and F } — 0, respectively. These forces are attractive for 
all interplate distances. 

Now we turn to the discussion of the general formula (|3.24|) for the vacuum forces 
in the limiting cases corresponding to small and large interplate distances. For this it is 
convenient to introduce a new integration variable v = y/t 2 — u 2 — m 2 as above and to 
pass to polar coordinates (r, 6) in the (u, v ) plane. For small distances the main contri¬ 
bution to the r-integral comes from large values of r. From the asymptotic behavior of 
the functions Fj(u) for large values of the argument described in section |2j. the functions 
Fj{u) in the coefficient of e 2at may be omitted and, to leading order, the corresponding 
forces coincide with the vacuum forces in the case of Neumann boundary conditions 
and are attractive. For large distance between the plates the main contribution to 
the integral results from small values of r and two subcases should be distinguished. 
For a massless scalar field or for a massive field with kernel functions fj(u) ~ u - ^', 
D/2 — 1 < j3j < D — 1, u — » oo, the terms with t in the coefficient of e 2at may be omitted 
and the vacuum forces coincide with those for the Dirichlet case and are attractive. In 
the second subcase, corresponding to the kernel functions fj(u) ~ o(u 1_£) ), u — > oo, to 


leading order the vacuum forces are given by expressions which are obtained from (|3.24|) 


by the replacements t ± F 3 (u) —> m ± Fj 0 in the coefficient of e 2at , where the constants 


Fj 0 are defined by formulae similar to (|2.9| ) with the replacements f(u) —> fj(u). 

We have evaluated numerically the vacuum forces acting on the plates in the case 
of the kernel functions 

./)(,•) - v (3.26) 

The corresponding Fourier transforms Fj(k\\) are given by formulae obtained from Eq. 


(|2.24|) by the replacements F\ — > f[ ] \ rj —> r]j. The parameters F^ are dehned by 


the formulae obtained from the corresponding expression for F\ in the paragraph after 
formula (|2.24|) by the replacements rj —>■ rjj, f 0 —> f 0 j. In figure 3 we have plotted the 
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vacuum pressure on the plate as a function on rja for a massless scalar in I) = 3 in the 
case of the kernel functions (|3.26|) with r]i = rj 2 = rj and F = 10. The numbers near 
the curves are the values of the parameter F\ 1 ^. For the values F^ < —1.08 the vacuum 
pressure is negative for all interplate distances and the corresponding vacuum forces 
are attractive. For the values F^ > —1.08 there are two values of the distance between 
the plates for which the vacuum forces vanish. These values correspond to equilibrium 
positions of the plates. For the values of the distance in the region between these 
positions the vacuum forces acting on plates are repulsive. Thus, the left equilibrium 
position is unstable and the right one is locally stable. 
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2 3 4 5 6 7 

Figure 3: Vacuum pressure on the plate, p/rj D+1 , as a function of ya for a massless scalar 
in D = 3 in the case of the kernel functions ( |3.2G| ) with r/i = rj 2 and = 10. The numbers 
near the curves are the values of the parameter i 7 ^ 1 ' defined in the paragraph after formula 

®)- 


19 












4. Concluding remarks 


In this paper we have investigated the positive-frequency Wightman function, the VEVs 
of the field square and the energy-momentum tensor for a scalar field with non-local 
boundary conditions on a single and two parallel plates in Minkowski spacetime. For 
the case of a single plate geometry we have considered the boundary condition of the 
form (|2.2| ), where the kernel function f(x n) describes the properties of the boundary. 
The VEVs of the physical quantities bilinear in the field are determined by the Fourier 
transform F(kn) of this function. By evaluating the corresponding mode-sum, we have 
presented the Wightman function as the sum of a Minkowskian part without boundary 
and boundary induced parts, formula (|2.10| ). The last term in this formula corresponds 
to the contribution of bound states. The boundary induced part in the VEV of the field 
square is obtained in the coincidence limit of the arguments of the Wightman function 
and is given by formula ( |2.14| ). The VEV of the energy-momentum tensor is obtained 
by acting with the corresponding second-order differential operator on the Wightman 
function and taking the coincidence limit. This VEV is determined by formula (|2.17|) . 
The vacuum stress in the direction orthogonal to the plate vanishes, and hence the 
corresponding vacuum force is zero. 

Unlike the case of local boundary conditions, in the non-local case the vacuum 
energy-momentum tensor does not vanish for a conformally coupled massless scalar 
field. Another important difference is that for non-local boundary conditions, in gen¬ 
eral, (T 0 °) sub ^ (Tp snh , n = 2,..., D. 

As in the case of local boundary conditions, the energy density and the vacuum 
stresses diverge on the surface of the plate. For a nonconformally coupled scalar field 
the leading term in the corresponding asymptotic expansion is the same as that for 
Neumann boundary condition. As an example, in figures 1 and 2 we have plotted the 
vacuum energy density and the vacuum stress as functions of the distance from the 
plate for a choice of kernel function / given by (|2.23| ) . 

In section [3] we have considered the geometry of two parallel plates with boundary 
conditions (Q. For the region between the plates the corresponding eigenvalues are 
solutions of equation (|3.5|) , where the coefficients Cj are determined by the Fourier 
transforms of the kernel functions fj(x\\) in the boundary conditions. The evaluation 
of the corresponding Wightman function is based on a variant of the generalized Abel- 
Plana summation formula, Eq. ( |3.12|) . The application of this formula allowed us to 
extract from the VEVs the parts resulting from the single plate and to present the part 
induced from the second plate in terms of integrals exponentially convergent for points 
away from the boundary. The Wightman function is presented in the form (|3.17|) . 
The VEVs of the field square and the energy-momentum tensor are obtained from the 
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Wightman function and are determined by formulae (|3.19|) , (|3.20| ). The vacnnm stress 
in the direction orthogonal to the plates is uniform. This stress determines the vacuum 
forces acting on the plates. The corresponding effective pressure is given by formula 

For small and large distances between the plates the vacuum forces are attractive. 
For intermediate distances the nature of the vacuum forces depends on the functions in 
the boundary conditions Q3.1p. For the example (|3.26|) we have shown that, depending 


on the parameters, the forces acting on the plates can be repulsive for intermediate 
distances (see figure 3). In this case it is possible to have a locally stable equilibrium 
value of the interplate distance stabilized by the vacuum forces. 

As far as we know, previous work on non-local boundary conditions for the Casimir 
effect had instead considered spectral boundary conditions for spinor fields in spheri¬ 
cally symmetric cavities [17], and hence our results with the non-local boundary con¬ 
ditions (2.2) and (3.1) for scalar fields are entirely original. From the point of view 
of bosonic string theory, a non-local Casimir effect is studied in |L8[, but of a com¬ 


pletely different nature as compared to our work, since one there deals with a non-local 
Lagrangian. 


Acknowledgments 

The work of A. Saharian has been supported by the INFN, by ANSEF Grant No. 05- 
PS-hepth-89-70, and in part by the Armenian Ministry of Education and Science, Grant 
No. 0124. The work of G. Esposito has been partially supported by PRIN SINTESI. 

References 

[1] Esposito G 1999 Class. Quantum Grav. 16 1113 

[2] Esposito G 1999 Class. Quantum Grav. 16 3999 

[3] D’Eath P D and Esposito G 1991 Phys. Rev. D 44 1713 

[4] Vassilevich D 2001 JHEP 0103 023 

[5] Grubb G 2003 Commun. Math. Phys. 240 243 

[6] Schroder M 1989 Rep. Math. Phys. 27 259 

[7] Esposito G 2000 Int. J. Mod. Phys. A 15 4539 

[8] Romeo A and Saharian A A 2002 J. Phys. A 35 1297 


21 


















[9] Saharian A A 2004 Phys. Rev. D 69 085005 

[10] Birrell N D and Davies P C W Quantum fields in Curved Space (Cambridge 
University Press, Cambridge, 1982) 

[11] Saharian A A 2006 Phys. Rev. D 73 044012 

[12] Deutsch D and Candelas P 1979 Phys. Rev. D 20 3063; Kennedy G, Critchley R and 
Dowker J S 1980 Ann. Phys. (N.Y.) 125 346; Milton K A 2004 J. Phys. A 37 R209 

[13] Ford L H and Svaiter N F 1998 Phys. Rev. D 58 065007 

[14] Graham N, Jaffe R L, Khemani V, Quandt M, Scandurra M and Weigel H 2002 Nucl. 
Phys. B 645 49; Graham N, Jaffe R L, Khemani V, Quandt M, Scandurra M and 
Weigel H 2003 Phys. Lett. B 572 196; Graham N and Olum K D 2003 Phys. Rev. D 
67 085014 

[15] Saharian A A 2000 The generalized Abel-Plana formula. Applications to Bessel 
functions and Casimir effect, Report No. IC/2000/14 [hep-th/0002239] 

[16] Saharian A A 2005 Bulk Casimir densities and vacuum interaction forces in higher 
dimensional brane models [hep-th/0508185]. 

[17] Cognola G, Elizalde E and Kirsten K 2001 J. Phys. A 34 7311 

[18] Elizalde E and Odintsov S D 1995 Class. Quantum Grav. 12 2881 


- 22 - 



